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ABSTRACT 

Narrow-Line Seyfert 1 (NLSl) galaxies show extreme properties with respect to the 
other Seyfert galaxies. Indeed, they are thought to be accreting at Eddington rates 
and to possess low mass black holes. Therefore, they may represent a key class of 
objects for understanding the co-evolution of black holes and their host galaxies. We 
propose that NLSls represent a class of AGN in which the black hole growth is, and 
has always been, dominated by secular evolution. Firstly, by looking at the NLSl host 
galaxy properties in the literature, we show that the evolution of NLSls is presently 
driven by secular processes, much more so than for Broad-Line Seyfert Is (BLSls). 
Secondly, we study the bulges of NLSl and BLSl galaxies. Our results demonstrate 
that NLSl host bulges are pseudo-bulges and are statistically different from BLSl 
bulges. This difference points to the particular importance of secular processes in the 
past evolution of their hosts. We build on this result to understand the implications 
on their evolution and the duration of their duty cycle. We show that NLSls are not 
necessarily in a special phase of black hole growth and that several Gyr are required 
for their black hole masses to become similar to BLSls. Finally, in the light of our 
results, we discuss the location of NLSl galaxies on the Mbh-ct plane and speculate 
about the connection between the NLSl galaxy properties and their black hole spin. 
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1 INTRODUCTION 

Since their discovery, Narrow-Line Seyfert 1 (NLSl) galax- 
ies have always been recognized as particular objects hold- 
ing important clues on the driving mechanisms of nuclear 
activity. First identified as objects similar to Seyfert Is with 
narrower Balmer lines, they were soon recognized as hav- 
ing exceptional spectral propertie s, both in their emission 
lines and in their continuum (see iKomossal |2008| . for a re- 
view). In fact, as shown through principal components anal- 
ysis (e.g. iBorosonlbOO^ I. NLSl galaxies are mostly clustered 
at one extreme end of the AGN parameter space (strongest 
FeII/[OIII] emission and lowest luminosity). Likewise, they 
are thought to possess small mass black holes and to have 
high Eddington accretion rates. In this sense, NLSls could 
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represent key objects in understanding the AGN phenomena 
and the co-evolution of massive black holes and their host 
galaxies. 

While the main defining criteria of NLSls with respect 
to BLSls is the empirical threshold at FWHM(H,3) ~ 2000 
km s~^, the properties of NLSls have been extensively stud- 
ied across many wavelength ranges. Trends and correla- 
tions have been identified using first s mall samples and later 
corroborated by larger su r veys (e.g. Veron-Cetty fc VeronI 
I2OOII : IWilhams et al.l l2002l : IZIiou et al.ll2006l '). Many scenar- 
ios have been considered to explain these propertie s, in par- 
ticula r their high accretion rates (L/I/Edd — 1, e.g. iBorosorJ 
l2002i : lGrupelli2004l . and reference therein) and low black hole 
masses (typically of order lO'^Mry^ , e.g. iBoller et all Il996l : 
IZhou et al.ll200d : iRvan et al.ll2007f ), but also, e.g., outflows, 
winds and density effects, high metallicity, p articular broad 
line region thicknesses and densities, etc. (see lKomossall2008l , 



2 G. Orban de Xivry et al. 



and references therein). While these scenarios can elucidate 
the nuclear properties of NLSls, they hardly explain the 
origin of the fundamental differences between NLSls and 
BLSls and, in particular, that NLSls appear to be more 
than just Seyfert Is with narrow lines. A few key questions 
could be formulated as follows: which particular mechanisms 
would lead to the Eddington accretion rates commonly seen 
in NLSls but observed less often in BLSls? What causes the 
difference in black hole growth of NLSls and BLSls that re- 
sults in low mass black holes in the former case? Could dif- 
fering host galaxy evolution explain the differences between 
NLSl and BLSl galaxies? 

In this paper, rather than studying the active nuclei, 
we investigate the host galaxies of NLSls and contrast their 
properties to those of BLSls, pursuing the hypothesis that 
different host galaxy evolution could explain the differences 
between NLSls and BLSls. In particular, we explore the 
relative role of secular processes in the evolution of NLSl 
and BLSl galaxies. Reviewing the literature on the mor- 
phology and the star formation in NLSl and BLSl hosts, 
we emphasize, in section [2j the present-day differences in 
their respective host galaxies. Afterwards, in section |3l we 
turn to the past evolution of NLSl and BLSl hosts. We 
perform a bulge-disk decomposition of samples of NLSl and 
BLSl galaxies and look at th eir bulge properties. Using pre- 
viously established c riteria /Kormo ndv fc Kennicutu l2004l : 
iFisher fc Drorvll2008l : [Cadotti 200a ). we are able to distin- 
guish pseudo- from classical bulges. This enables us to de- 
termine the main processes that have driven the evolution 
of the NLSl and BLSl hosts. We analyse the differences of 
NLSl and BLSl host bulge property distributions, conclud- 
ing that NLSl galaxies contain pseudo-bulges and, hence, 
have always been dominated by secular evolution. Finally, 
in section [3] and section [S] we investigate the cosmological 
context of the NLSl host phenomenon driven by such an 
evolutionary mode. We then note the link between secular 
evolution and rapidly spinning black holes, and speculate on 
the location of NLSl galaxies on the Mbh — cr relation. We 
conclude by summarizing our picture of the NLSl galaxy 
phenomenon and present ways to further test our proposi- 
tion. 

When calculating distances and look-back times, we 
assume a general relativistic Friedmann-Robertson- Walker 
(FRW) cosmology with matter-density parameter f2m = 0.3, 
vacuum energy- density parameter Q.a = 0.7 and Hubble 
parameter Hq = 70 km s~^ Mpc~^. 



2.1 Morphological properties 

2.1.1 Large-scale bars 

The morphology o f NLSl host galaxies has been s t udied 
in several paper s llCrenshaw et al.l l2003l : iDeo et al.l l2006l : 
lOhta et al.l l2007l : iRvan et al.ll2007^ . The main results can 
be summarized as follows: 

NLSl host galaxies are likely to he strongly barred (much 
more than BLSl ones) and their nuclear dust morphology is 
likely to be a grand-design spiral. 

The bar frequen cy among NLSl and BL Sl host galaxies 
has be en studied bv lCrenshaw et al] ||200^ ) and Ohta ct al] 
|2003). The first paper reports a visual study based on an 
HST survey of 91 Seyfert galaxies (13 NLSls and 78 BLSls) 
at 2 < 0.035 (Mal kan et al] 1 19981 ) and 6 additional NLSls 
at z < 0.084 (HST archives in the IVeron-Cettv et allboOlal 
sample). The results are striking: 91% of the sample is classi- 
fied as spiral galaxies among which 65% of NLSls have bars, 
and 25% of BLSls have bars. More particularly, the authors 
look at the fraction of barred spiral galaxies in their sample 
as a function of the full width at half maximum (FWHM) 
of the broad component of the H,3 emission line. As pre- 
sented in Fig. [T] they obtain a clear difference between the 
two regimes of NLSls and BLSls. We note that in this fig- 
ure, since FWHM measurements are not available for every 
BLSl, we have made a single bin for th e BLSl class. 

The results from lOhta et al.l (|2007| ) are more conserva- 
tive. They use an heterogeneous set of data of NLSl galaxies 
and look at the morphology and the possible trends with the 
NLSl properties. They perform a visual and a quantitative 
classification based on ellipse fitting of isophotes. While they 
confirm the high bar fraction among (spiral) NLSl hosts, 
they do not observe a clear trend with the FWHM. Never- 
theless, if we consider only the fraction of spirals with strong 
bars (SB but not SAB), the trend does appear clearly using 
their visual classification (as represented in Fig.[T|. Turning 
to their quantitativ e classification, we note that one of the 
criteria lOhta et al.l use to identify the bars is an ellipticity 
Ebar ^ 0.25, whcrc Ebar = niax(egaiaxy ) • Howcver a common 
practice to identify strong bars is to use ebar ^ 0.45 (e.g. 
[Shi osman et al.ll2000l ). Applying this latter criterion on their 
sample by examining the ellipse fit plots, we obtain, for the 
respective bins in Fig. [T] (i.e. 500-1000, 1000-1500, 1500- 
2000 km s"^), bar fractions (i.e. 89%, 46% and 57%) similar 
to their visual classification (i.e. 100%, 50% and 57%, see 
Fig.HI. 



2 SECULAR EVOLUTION IN NLSl HOST 
GALAXIES 

In this section, we review the literature that has been pub- 
lished concerning the differences in the host galaxy prop- 
erties of NLSls and BLSls. We focus in particular on the 
morphology and the star formation rate, and emphasize the 
role of present secular processes in distinguishing between 
these two classes of type 1 AGN. 



2.1.2 Circumnuclear morphology 

In a study of nuclear dust morphology in a matched-paired 
sample of a ctive/inact ive galaxies and barred/unbarred 
galaxies. Mar tini et al.l (2003) show that grand-design nu- 
clear dust spirals are only found in galaxies with a large-scale 
bar. However, while not finding any universal nuclear mor- 
phology in active galaxies, they do find similar features in 
the circumnuclear environ ments of bo th active and inactive 
galaxies. In another study. [Peo et alj (|2006) investigate the 
nuclear dust morphology in NLSl and BLSl host galaxies 
based only on the HST survey conducted by iMalkan et al.l 
(1998). Their study also shows that the grand-design nu- 
clear dust spirals are largely present in barred galaxies. They 
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Eddington ratio and stronger star formation activity com- 
pared to their broad-line counterparts. 



2.3 



Secular processes all the way to the SMBH 



As discussed above, the current morphology of NLSl host 
galaxies is distinguishable from other Seyfert galaxies. In- 
deed, in contrast to BLSls, NLSl galaxies are likely to be 
strongly barred and to show more intense central star for- 
mation. This is in line with the fact t hat bars are known to 
drive gas into the central kiloparsecs (jSakamoto et al.lll999l : 
ISheth et al]|2005h . and that nuclear star formation is en- 
hanced in barred galaxies ||Ho et al.lll997l ). 

While no universal fueling mechanism for low- 
luminosity AGNs seems to operate in galactic nuclei 
(jMartini et al.ll2003l ). the NLSl host morphology typically 
exhibits a circumnuclear grand-design spiral. This appears 
to be linked to the p resence of strong bars l|Martini et alj 
l2003l : lDeo et al.|[200ff ). and indeed bars are able to drive cir- 
cumnuclear spiral structures ( Maciejewskj ,2004a ibh ■ Hence, 
NLSl galajcies show uninterrupted asymmetries able to drive 
the gas inwards from a few kpc to a few tens of pc. The par- 
ticular strength of secular processes in NLSls could therefore 
account for the high central star formation and presumably 
to the large Eddington rates observed in NLSls. 



FWHM [kms-i] 



Figure 1. Histograms presenting the fraction of Seyfert 1 spirals 
a) with bars and b) with nuclear grand-design spirals as a func- 
tion of the FWH M of their broad emissi o n line, Hfi. Drawn from 
data reported in ICrenshaw et al et al ] l|2007l 'l and 

iDeo et al.l ||2006| '). 



classify the nuclear structures and find that (1) the nu- 
clear dust morphologies in NLSl/BLSl, Barred/Unbarred 
are mainly nuclear dust spirals, (2) in the "nuclear dust 
spiral" class, NLSls are more likely to have grand-design 
spirals than BLSls. Fig. [T] also shows that this fraction of 
grand-design spirals in NLSls follows the same trend with 
the FWHM(H,3) as the bar fraction. 

As we exp ect strong bars to d rive a circumnuclear spi- 
ral str ucture (iMacieiewskil 2004a ^ and to d rive gas in- 
wards (jSakamoto et al.lll999l : ISheth et al.ll2005l l. we expect 
the presence of such asymmetries in host galaxies to result 
in an enhanced star formation in the central kiloparsecs. 



2.2 Star formation 

In a recent paper, ISani et ahl (|2010l ^ study the link between 
star formation in the central kpc vs. FWHM(H;3) in NLSl 
and BLSl host galaxies. After discussing carefully possible 
luminosity and distance effects, they conclude that NLSls 
are associated with more intense star formation than BLSls 
(with on average a star formation to AGN ratio > 2 times 
larger in NLSls). More generally, they find that type 1 AGNs 
with narrower broad emission lines reside in hosts containing 
more intense star-forming regions. 

Finally, they find a connection between high Eddington 
ratio and high star formation rates concluding that NLSls 
are characterized by smaller black hole (BH) mass, larger 



3 BULGES OF NLSl HOST GALAXIES 

Since strong secular evolution is currently occurring in 
NLSls, it is important to ask whether or not secular pro- 
cesses have shaped the NLSl host galaxies by dominating 
their past evolution and hence influencing their black hole 
growth. 

We address this issue by examining the bulges of NLSl 
host galaxies and comparing them to those of BLSls, since 
one can expect to observe evolutionary dependent bulge 
characteristics. Specifically, an evolution driven mainly by 
galaxy mergers will result in different bulge properties than 
if the evolution is mainly driven by internal secular evolu- 
tion. 

In this section, we compare NLSl and BLSl galax- 
ies by performing a photometric bulge-disk decomposition 
of homogeneous samples of NLSls and BLSls. This com- 
parison is put in perspective with previous studies on the 
distinction b etween pseudo- and classical bulges in inac- 
tive galaxies |Kormendv fc Kennicuttll2004l : iFisher &: Drorvl 
120081 : iGadottl I2009I '). and also on the bulge-disk decomposi- 
tion of two galax: y samples composed exclusively of NLSls 
[Mathur et al.ll20Tll ). 



tion of two galax:y i 
(|Rvan et al.ll2007i : 



3.1 Pseudo-bulges and secular evolution in disk 
galaxies 

iKormendv fc KennicuttI (|2004h review in detail the forma- 
tion of pseudo-bulges by secular processes. As dense central 
components of galaxies, pseudo-bulges differ from classical 
bulges in that they were made slowly by disks out of disk 
material while classical bulges are "merger-built" bulges. 
Therefore, pseudo-bulges are formed by internal secular pro- 
cesses such as bar instabilities, spiral structures, etc. as op- 
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posed to galaxy mergers or external secular evolution (minor 
mergers, prolonged gas infall, etc.). 

As pseudo-bulges retain memory of their disky origin, 
it is possible to disent a ngle them from classical bulges. 
iKormendv &: KennicuttI l|2004 ) pointed out the Sersic in- 
dex as one way to identify them. Indeed, since a pseudo- 
bulge forms from gas accreting from the disk, it has a sur- 
face brightness profile similar to that of the outer disk and 
ther efore would have a low Sersic index n;, ~ 1 — 2. 

iFisher fc Drorvl (|2008l ') have studied in great detail the 
structure of classical bulges and pseudo-bulges using high- 
resolution data (77 inactive galaxies with data in the HST 
archive and z <0.01). They use morphological signatures to 
first visually classify the bulges of nearby galaxies as pseudo 
or classical. They then perform a bulge-disk decomposition 
and study in particular the distribution of Sersic indices. 
They find that, statistically, pseudo-bulges have Sersic in- 
dices rij, < 2 while classical bulges have Uf, > 2. This result 
shows that the Sersic index is a good statistical tool to test 
if a class of objects has cla ssical or pseudo-bulges. 

Finally, fCadottl (|2009l ) also study pseudo- and classical 
bulges using a large, low resolution, SDSS sample of galax- 
ies (~1000 inactive galaxies with 0.02 ^ 2: ^ 0.07). Ife uses 
the position in the < fi^ > —re plane, also called the Ko- 
rmendy relatiorQ, to study the bulge properties and iden- 
tifies pseudo-bulges as being fainter in surface brightness 
for a given half-light radius (much fainter than predicted 
by the correlation fit to elliptical galax;ies). Where iGadottil 
l|2009h clearly sees independent groups in his i-band den- 
sity plot of the Kormendy r elation (< jie > —r^ relation), 
iFisher fc Drorvl (|2008l . I2OIOI ) only find that pseudo-bulges 
scatter around the photometric projections of the funda- 
mental plane. 

Neverthe l ess, all these authors (|Fisher fc Drorvl l2008l : 
lGadotti|[2009l : IFisher fc Dror^l2010l ') agree that most of the 
pseudo-bulges have a low Sersic index rii, < 2 and that they 
tend to be less prominent than classical bulges (in particular, 
they tend to have a low bulge-to-total light ratio). 

Based on these considerations for bulge classification in 
inactive galaxies, we will use the Sersic index to identify the 
prevailing bulge type in the NLSl and BLSl populations as 
pseudo- or classical. We will then study the prominence of 
NLSl and BLSl bulges. 



3.2 Bulge/disk decomposition 

Building on the work of IFisher fc Drorvl (|2008l ) for inactive 
galaxies, we select archive HST images of Seyfert galaxies 
to study the bulges of active galax;ies, in order to assess 
whether the bulge characteristics of NLSls and BLSls might 
explain, by their evolutionary implications, the distinctions 
between these two classes of AGN. Crucially, by performing 
the bulge-disk decomposition for samples of both NLSls and 
BLSls, we minimize the impact of any systematic errors that 
our fitting procedure might generate. 

We sele c t NL Sl and BLSl galaxies from the 
iMalkan et all (| 19981 ) HST imaging survey of nearby AGNs. 



^ relation between the mean effective surface brightness within 
the effective radius < fie > and the half-light radius re, which is 
a projection of the photometric fundamental plane. 



This survey contains a uniform sample of 91 Seyfert 1 galax- 
ies at 2: ^ 0.035 observed with the Wide Field Planetary 
Camera 2 ( WFPC2) through the F606 W fiher. T his sample, 
also used bv lCrenshaw et ahl (|2003. ) and lDeo et al . (2006) in 
their morphological studies, contains 11 NLSl galaxies and 
80 BLSl galaxie s. The 11 NLSl galaxies ar e genuine NLSls 
as identified by IVeron- Cettv et all (|2001bl ) on the basis of 
their optical spectra (broad component of Hj3 < 2000km 
s~^ a nd strong FeH emissio n ) and are listed in the cata- 
log of|v eron- Cettv fc Veroiil ((2OI0I ). We therefore initially 
select all the 11 NLSl galaxies available, as well as 21 of 
the 80 BLSls. The BLSl sample selection is made in a 
way to roughly matc h the ^25% fraction o f such hosts that 
are strongly barred l|Crenshaw et al.ll2003l ). Their individ- 
ual Seyfert classifications are reported i n Tab.[T] These are 
taken from lVeron-Cettv fc VeronI (I2OI0I) with 2 exceptions. 
IC1816 is classified by [Nlarqucz et al.l (|2004l ') as a type 1. 
Because their spectrum shows clear evidence of broad emis- 
sion, we have adopted this classification. For NGC5252, 
we have followed the classi fication as a type 1.9 given in 
lOsterbrock fc Martell (|l993h . because the presence of broad 
H a, with a measured FWHM of ~2500 km s"'^, is confirmed 
bv lAcosta-Puhdo et al.1 (| 19961 '). We iiote th at the broad Ha 
is very clear in polarized light (|Tranll20"lol ). and shows dra- 
matic variations over a period of several years. This may 
be why there is some uncertainty about its classification as 
a type 1 or 2. Following this initial selection, we refine our 
sample based on the limited field of view (FoV) of WFPC2 
by rejecting objects with z < 0.010 or scales ^0.23 kpc/" in 
order to obtain a reasonable minimum FoV of >8x8kpc on 
each object. This criterion ensures that the disk of the host 
can be fitted properly, while the redshift limit of the orig- 
inal source catalogue ensures that the bulge is sufficiently 
well resolved. The two final samples for which we performed 
the bulge-disk decomposition are composed of 10 NLSls and 
19 BLSls, and are given in Tab. \T\ Finally, we check that 
the mean redshift is not strongly biased with respect to the 
NLSl sample (< 2:nlsi >= 0.024, < 2:blsi >= 0.027), and 
that no particular circumnuclear morphology has been se- 
lected. 

To perform the bulge-disk decomposition of these 29 
galaxies, we use the two-dim ensional profile fitting algorithm 
GALFn0 dPeng et al.ll2010l ). For each galaxy we iteratively 
fit three components: a Gaussian profile, a Sersic profile and 
an exponential profile (see also Appendix [X]) . These com- 
ponents are aimed at modelling respectively the nuclei, the 
bulges and the disks of our galaxies. While these HST data of 
nearby AGN have small pixel scales (see column 4 in Tab.[T]) 
enabling us to better constrain the central regions of the 
galaxies, the FoV of WFPC2 (~ 35"x35") is smaU relative 
to the full extent of the galaxies. This makes it hard to con- 
strain the sky level. This issue has already been addressed 
to some extent by the refinement of our sample selection to 
objects with z ^ 0.010. We cover it further in Appendix [Bl 
when we discuss our treatment of the possible coupling be- 
tween the background level and the exponential profile. Fi- 
nally, we analyze the robustness of our fit by studying the 
effect of saturated regions in the images and the dependence 
of the fits to the PSF. We detail our iterative fit procedure 
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Table 1. HST sample of NLSls and BLSls. Columns: (1) object name; (2)-(4) J2000 coordinates and redshift from Nasa/IPAC Ex- 
tragalactic database (NED); (5) luminosity distance in Mpc for an Hq = 70km s~^ Mp c~^, Vlrn = 0-3 and = 0.7 cosmology; (6) 
respective scale, the WF PC2 pixel scale is 0.0456"; ( 7) morphological classification (from lMalkan et al.lll99^ . MGT); (8) Seyfert classi- 
fication according to the IVeron- Getty fc VeronI ||2010| ') catalog except for IC1816 and NGC5252, see also the main text. SI classification 
designates a type 1 A GN with unspecified sub-ty pe; (9) FWHM of the broad component of (or in a f ew ca ses, Ha), the NLSl 
measurements are from lv^ron-Cettv et al ] l|2001bl ). the BLSl measurements are taken from [Crenshaw et al ] l|2003h . see also references 
therein. 



Object Name R.A. Dec. z D px sc. Morpho. AGN FWHM 

(J2000) (J2000) (Mpc) (kpc/") (MGT) type (km s-i) 



NLSl SAMPLE 


KUG1136 


11 


39 


13.9 


-1-33 


55 


51 


0.032 


131.84 


0, 


.64 


SBO 


Sin 


1145 


MRK0042 


11 


53 


41.8 


-1-46 


12 


43 


0.024 


99.83 


0, 


.48 


SBa 


Sin 


865 


MRK0335^ 


00 


06 


19.5 


-1-20 


12 


11 


0.025 


109.12 


0, 


.50 


7 


Sin 


1350 


MRK0359"~ 


01 


27 


32.5 


+ 19 


10 


44 


0.017 


71.31 


0, 


.35 


SBb/c 


Sin 


900 


MRK0382 


07 


55 


25.3 


+m 


11 


10 


0.034 


139.74 


0, 


.68 


SBa 


Sin 


1280 


MRK0493 


15 


59 


09.6 


-1-35 


01 


47 


0.031 


127.87 


0, 


.62 


S(B)a 


Sin 


740 


MRK0766 


12 


18 


26.5 


4-29 


48 


47 


0.012 


50.65 


0, 


.25 


SBc 


Sin 


1630 


MRK0896 


20 


46 


20.8 


-02 


48 


45 


0.027 


111.91 


0, 


,54 


Sc 


Sin 


1135 


MRK1044 


02 


30 


05.5 


-08 


59 


53 


0.016 


67.20 


0, 


,33 


Sa 


Sin 


1010 


NGG4748 


12 


52 


12.4 


-13 


24 


53 


0.014 


58.94 


0, 


,29 


Sa 


Sin 


1565 


BLSl SAMPLE 


ES0438G9£ 


11 


10 


48.0 


-28 


30 


04 


0.024 


99.83 


0, 


.48 


SBc/d 


SI 


5000 


F1146 ~ 


08 


38 


30.8 


-35 


59 


33 


0.032 


131.84 


0, 


,64 


Sb 


SI 


4300 


IC1816 


02 


31 


51.0 


-36 


40 


19 


0.017 


71.31 


0, 


,35 


SBa/b 


SI 




Mrk0279 


13 


53 


03.4 


-1-69 


18 


30 


0.031 


123.89 


0, 


,60 


Sa 


Sl.O 


6860 


Mrk0290 


15 


35 


52.3 


+57 


54 


09 


0.029 


123.89 


0, 


,60 


E 


SI. 5 


2550 


Mrk0352 


00 


59 


53.3 


-h31 


49 


37 


0.015 


63.08 


0, 


,31 


E 


Sl.O 


3800 


Mrk0423 


11 


26 


48.5 


-f35 


15 


03 


0.032 


131.84 


0, 


,64 


Sb 


SI. 8 


9000 


Mrk0530 


23 


18 


56.6 


-fOO 


14 


38 


0.029 


119.91 


0, 


,58 


Sa 


SI. 5 


6560 


Mrk0595 


02 


41 


34.9 


4-07 


11 


14 


0.028 


111.91 


0, 


,54 


Sa 


SI. 5 


2360 


Mrk0609 


03 


25 


25.3 


-06 


08 


38 


0.032 


139.74 


0, 


,68 


Sa/b 


SI. 8 




Mrk0704 


09 


18 


26.0 


+ IQ 


18 


19 


0.029 


119.91 


0, 


,58 


SBa 


SI. 2 


5500 


Mrk0871 


16 


08 


36.4 


+ 12 


19 


51 


0.034 


139.74 


0, 


,68 


Sb 


SI. 5 


3690 


Mrk0885 


16 


29 


48.2 


-1-67 


22 


42 


0.026 


103.87 


0, 


,50 


SBb 


Sl.O 




Mrkll26 


23 


00 


47.8 


-12 


55 


07 


0.010 


46.48 


0, 


,23 


Sb 


SI. 5 




Mrkl400 


02 


20 


13.7 


4-08 


12 


20 


0.029 


119.91 


0, 


,58 


Sa 


Sl.O 




NGG5252 


13 


38 


15.9 


4-04 


32 


33 


0.022 


95.79 


0, 


,46 


SO 


SI. 9 


2500 


NGG5940 


15 


31 


18.1 


4-07 


27 


28 


0.033 


139.74 


0, 


,68 


SBc 


Sl.O 


5240 


NGG6212 


16 


43 


23.1 


4-39 


48 


23 


0.030 


123.89 


0, 


,60 


Sb 


SI 


6050 


IISZIO 


13 


13 


05.8 


-11 


07 


42 


0.034 


139.74 


0, 


,68 


7 


SI. 5 


3760 



From the bulge disk decomposition, MRK335 seems better described by a unique high Sersic profile. This 
object is therefore rejected from our bulge analysis. 

^ From the bulge disk decomposition, ES0438G9 seems a bulgeless galaxy. This object is therefore rejected 
from our bulge analysis. 



in Appendix |B] our treatment of additional structures such 
as bars, rings and spirals, and the particular attention given 
to the background level. We give also eight examples of our 
fits in Fig. [B2] and Fig. [H 

The relevant results of our fits are given in Tab. [2] 
Among them, the reduced chi-square value from the fit, 
first indicator of its quality. Our mean values are 1.1 for 
NLSls and 1.06 for BLSls, reflecting the overall acceptabil- 
ity of the fits. 

During this process, we have found two objects for 
which no acceptable bulge disk decomposition could be per- 
formed. Specifically, the morphology of MRK335 is mainly 
point-like and is better described by a single high Sersic in- 
dex profile, while ES0438G9 seems to be a bulgeless galaxy 



also consistent with its morphological type. We have there- 
fore excluded these two objects from the bulge analysis dis- 
cussed in Section r3.3l which is performed on 9 NLSl and 18 
BLSl galaxies. 

3.3 Structural properties of NLSl and BLSl host 
galaxies 

3. 3. 1 The Sersic index rib in the bulges 

In Fig. (2] we compare the results o btained for the 9 NLSl 
galaxies to the distribution found bv lFisher fc Drorvl (|2008l ) 
for pseudo- and classical bulges. According to their results, 
the mean Sersic index of pseudo-bulges is 1.69 with only 
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Table 2. Results of the bulge disk decomposition for the NLSls and BLSls. Columns: (1) object name; (2) fitted components. The 
different components p, g, s, d, b stand respectively for PSF, Gaussian, Sersic, disk (exponential) and background (sky). Components 
are put in brackets if one or more of their parameters are kept fixed in the fit; (3) the FWHMg of the Gaussian component in kpc; (4) 
the bulge Sersic index; (5) i?;,, the effective radius of the bulge in kpc; (6) Ra, the scale radius of the disk in kpc; (7)- (8) axis ratio of 
the bulge and the disk components; (9)-(10) B/D and B/T, the bulge-to-disk and bulge-to-total luminosity ratios; (11) x5i the reduced 
of the obtained fit. See also Appendix [Al 



Object Name 


Comp. 


FWHMg 


rib 


Rb 




lb 


Id 


B/D 


B/T 


X.1J 






(kpc) 




(kpc) 


(kpc) 












NLSl SAMPLE 


KUG1136 


g+s+d-l-b 


0.07 


1.21 


0.95 


3.15 


0.76 


0.62 


0.20 


0.17 


1.16 


MRK0042 


g+s+d-l-b 


0.07 


1.27 


0.44 


3.12 


0.82 


0.58 


0.21 


0.18 


0.95 


MRK0335^ 


g+s+b 


0.18 


3.41 


1.27 




0.93 








1.20 


MRK0359~ 


g+s+d-(-b 


0.08 


1.41 


1.19 


7.45 


0.67 


0.83 


0.12 


0.11 


1.03 


MRK0382 


g-fs-fd-(-b 


0.12 


1.36 


0.52 


3.02 


0.94 


0.54 


0.46 


0.31 


1.49 


MRK0493 


g-l-s-fd-(-b 


0.08 


0.74 


0.34 


3.74 


0.97 


0.41 


0.17 


0.15 


1.02 


MRK0766 


g-fs-t-d-fb 


0.06 


1.88 


0.15 


1.21 


0.74 


0.44 


0.13 


0.11 


1.05 


MRK0896 


g-fs-fd-fb 


0.10 


2.06 


0.37 


2.82 


0.77 


0.71 


0.18 


0.15 


1.34 


MRK1044 


g-l-s-l-d-l-b 


0.07 


1.45 


0.20 


1.18 


0.91 


0.76 


0.44 


0.30 


1.05 


NGC4748 


g-l-s-|-d-(-b 


0.07 


1.93 


0.25 


1.89 


0.98 


0.68 


0.25 


0.20 


0.68 


BLSl SAMPLE 


ES0438G9^ 


g-l-d-l-b 


0.06 






1.76 




0.52 






1.39 


F1146 ~ 


g-l-s-l-d-l-b 


0.07 


3.74 


0.90 


2.09 


0.47 


0.65 


1.22 


0.55 


1.54 


IC1816 


g-|-s-|-[d]-|-b 


0.07 


1.98 


0.41 


[4.47] 


0.92 


0.67 


0.11 


0.10 


0.84 


Mrk0279 


p-|-s-|-d-|-[b] 




2.18 


2.48 


11.11 


0.58 


0.56 


0.57 


0.36 


0.97 


Mrk0290 


g-l-s-l-d-l-b 


0.10 


4.06 


0.47 


2.22 


0.86 


0.88 


0.90 


0.47 


1.11 


Mrk0352 


g-l-s-l-d-l-b 


0.05 


4.49 


0.90 


1.59 


0.98 


0.76 


0.79 


0.44 


0.84 


Mrk0423 


[s-Hd]-Hb 




2.13 


0.43 


1.72 


0.69 


0.68 


0.73 


0.42 


1.62 


Mrk0530 


[g]-l-s-|-d-|-b 


0.17 


2.4 


1.04 


4.04 


0.85 


0.65 


0.60 


0.38 


0.75 


Mrk0595 


g-f-s+[d]+b 


0.06 


3.47 


0.74 


1.86 


0.66 


[0.75] 


1.92 


0.66 


1.30 


Mrk0609 


g-fs-f[d]-f[b] 


0.15 


2.28 


1.48 


[2.53] 


0.79 


[0.95] 


1.60 


0.62 


2.52 


Mrk0704 


g-l-s-l-d-l-b 


0.16 


2.88 


1.19 


6.29 


0.70 


0.50 


0.75 


0.43 


1.25 


Mrk0871 


g-l-s-l-d-l-b 


0.12 


1.28 


0.61 


3.78 


0.52 


0.38 


0.13 


0.11 


0.91 


Mrk0885 


g-l-s-l-d-l-b 


0.07 


2.62 


2.70 


9.05 


0.74 


0.52 


0.32 


0.24 


0.83 


Mrkll26 


g-l-s-l-d-l-b 


0.04 


1.86 


0.27 


1.78 


0.86 


0.66 


0.17 


0.15 


0.47 


Mrkl400 


g-l-s-l-d-l-b 


0.12 


1.7 


0.55 


2.27 


0.61 


0.31 


0.32 


0.24 


0.98 


NGC5252 


g-l-s-l-d-l-b 


0.07 


3.9 


2.42 


3.20 


0.53 


0.44 


0.58 


0.37 


0.35 


NGC5940 


g-l-s-l-d-l-b 


0.10 


1.23 


0.30 


3.77 


0.89 


0.66 


0.06 


0.06 


0.99 


NGC6212 


g-l-s-l-d-l-b 


0.09 


1.52 


1.23 


2.88 


0.75 


0.70 


0.77 


0.43 


0.87 


IISZIO 


g-l-s-l-d-l-b 


0.18 


1.92 


1.56 


4.29 


0.91 


0.72 


0.70 


0.41 


0.57 



From the bulge disk decomposition, MRK335 seems better described by a unique high Sersic profile. This 
object is therefore rejected from our bulge analysis. 

^ From the bulge disk decomposition, ES0438G9 seems a bulgeless galaxy. This object is therefore rejected 
from our bulge analysis. 



~10% of them having an index greater than 2. For our NLSl 
sample, we find < Ub >~ 1.48 (and a standard deviation for 
the distribution of (t„ ~ 0.39) with none of them signifi- 
cantly exceeding a Sersic index of 2. 

The bottom panel of Fig. [2] compares the distribution 
of NLSl and BLSl host bulge Sersic indices from the anal- 
ysis of our samples. While the two distributions are clearly 
different, BLSl host bulges also tend to have lower Sersic 
indices than the lFisher fc Drorvl (|200^ classical bulges. We 
find < Ub >~ 2.54, (t„ ~ 0.97 for B LSl host bulges in 
contrast to < Ub >~ 3.49 obtained bv lFisher fc Drorvl for 
classical bulges. Thus, our results suggest that BLSls do 
not have "pure" classical bulges, but rather mixed bulges 
composed of pseudo- and classical components. 



It is approp r iate to mention here the work of 
iLaurikainen et alj (|2007h . They found that, for inactive 
galaxies, the mean bulge Sersic index is ~2.5 or less across 
the Hubble sequence. These results can also be interpreted 
as the existence of a large range of composite bulges between 
the two extreme "pseudo-" and "classical" bulge types. 
Putting this result into perspective with our bulge/disk de- 
composition suggests that the bulges of NLSl hosts are likely 
to be "pure" pseudo-bulges, while the bulges of BLSl hosts 
appear to be composite bulges and hence have Sersic indices 
distributed around Ub ~ 2.5. 

In order to test the significance of the difference be- 
tween the Sersic index distributions of NLSls and BLSls, 
we use the Kolmogorov-Smirnov test. As such, the cumula- 
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Figure 2. Histogram of bulge Sersic indices, rtj,. a) NLSl host 
bulges (9 objects) from our sample compared to pseudo-bulges 
(53) and classical bulges (26, we do not includ e their sample 
of elliptical galaxies) from iFisher fc Drorvl | |2008|) . b) NLSl host 
bulges compared to BLSls host bulges (18 objects). 



tive distributions of the Sersic index are presented in Fig. [S] 
The slopes of these distributions reflect the difference in the 
dispersions given above ~ 0.39 and 0.97 for the NLSl 
and BLSl samples respectively). They emphasize that the 
population of NLSl bulges appears to be like the popu- 
lation of pseudo-bulges, while the properties of the BLSl 
bulges are more widely distributed between pseudo- and 
classical bulges. While we cannot fully reject the null hy- 
pothesis of the Kolmogorov-Smirnov test, it yields a proba- 
bility of < 1.2% that the NLSls and the BLSls are drawn 
from the same parent distribution. This is a remarkable re- 
sult that again underlines the connection between black hole 

and bulge properties. 

Since a result from lWeinzirl et al.l l|2009l ) suggests a pos- 
sible link between bars and bulges, we verify that the rela- 
tive bar fractions in our samples have a minimum impact on 
our result of Fig. [2] Indeed, as described previously in sec- 
tion 13.21 our selection of NLSls and BLSls respects the rel- 
ative bar fraction observed, i.e. ~ 25% in BLSl and ~ 75% 
in NLSl galaxies. Therefore, we consider in Fig. |4] (right) 
the distribution of bulge Sersic indices of NLSl and BLSl 
galaxies with bars only. While we observe a general small 
shift towards lower Sersic indices (as expected from, e.g., 
IWeinzirl et~aLll2009l ), NLSl and BLSl are clearly distinct. 
Following this observation, we reasonably conclude that al- 
though the presence of a bar can be linked to the bulge Sersic 
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Figure 3. Cumulative distr ibution vs. the Sers ic index for 
pseudo- and classical bulges jFisher fc Drorvl |2008|) . and NLSl 
and BLSl host bulges, this paper. 



index, it does not imply that the bulge is a pseudo-bulge. 
Thus, our result that NLSls tend to possess pure pseudo- 
bulges in contrast to BLSls, is not a consequence of different 
bar fractions in the two populations. 

Finally, we consider the influence of the Seyfert type 
of the BLSls on our Sersic measurements. Indeed, changes 
from type 1 to intermediate type Seyferts can be attributed 
to changes in the ionizing radiation of the AGN (e.g. 
lGoodrichlll990D. or to variation in t he absorbing material 
('Maiolino & Rieke"l995'; 'Forstei"l999') , or as well to different 
inclination of the host galaxies (e.g. Maiolino & Rickc 199^; 
iRamos Almeida et al.l 120091 ). These effects can explain the 
large Balmer decrements in the optical spectra of intermedi- 
ate type galajcies but can also indirectely bias our Sersic in- 
dex measurements. Indeed, in the case of intermediate type 
Seyferts, a fainter AGN would ease the fit, but an excess of 
dust in the host galaxy and projection effects (nearly edge-on 
host spirals) can make fitting the host harder. For these rea- 
sons, we check whether there is any systematic effect on our 
fits by looking at the Sersic index versus the Seyfert type of 
the BLSls (given in Tab.[T|, the result is given Fig.\^(left). 
Since no trend of the Sersic index with the BLSl sub-class 
can be observed, we can be confident that our fits are not bi- 
ased by systematic effects related to the intermediate BLSl 
classifications. 

Similarly, NLSls are a sub-class of type 1 AGN, for 
which the Seyfert sub-classification is based solely on the 
properties of the AGN itself. But our result also shows that 
the bulge properties of the NLSl host galaxies do represent 
a distinct sub-class of the bulges of BLSl hosts. 



3.3.2 Bulge prominence and fundamental plane 
projections 

Building on our fit results, we study the bulge prominence 
by looking at the size distribution of the bulges, as well 
as the bulge-to-disk (B/D) and bulge-to-total (B/T) light 
ratios. The first objective here is to further explore the dif- 
ferences between NLSl and BLSl host bulges. But doing 
so also enables us to confirm the validity of our bulge-disk 
decomposition. 
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Figure 4. (Left) BLSl Seyfert type versus their respective Sersic 
index. Since there is no trend, the Sersic index does not seem 
affected by the Seyfert classification of our BLSls. (Right) His- 
togram of bulge Sersic indices, n;,, of NLSl (6 objects) and BLSl 
(4 objects) barred galaxies, see also Tab. [T] for the classification 
and Tab. [2] for n^. 



We compute the B/D and B/T light ratios using our 
fit parameters (see Equation IA3I and Equation IA4|) given 
in Tab. [51 and present their distribution in Fig. [S] (left). 
The median B/T of 0.39 in BLSls and 0.17 in NLSls indi- 
cates that NLSl galaxies have lower B/T light ratios than 
BLSls. We compare these distributions to <B/T>= 0.41 
for an average classical bulge a nd <B/T>= 0.16 f or an av- 
erage pseudo-bulge given by Fis her fc Drorvl l|2008l ). [Gadottil 
l|2009 l also finds similar values. The B/T ratio therefore pro- 
vides strong support for our conclusion that NLSl bulges are 
pseudo-bulges, while BLSl bulges are largely composite or 
classical. 

In the same Fig. [5] (right), we also plot the B/T 
light ratios versus the Hubble Type, the Sersic index, and 
the effective radius of the bulges. As expected, the mean 
B/T ratio tends to decrease with the Hub ble Type (e.g. 
iGraham fc Worleylliooi : [Masters et al.ll2010l '). and hence it 
appears that NLSl galaxies tend to be of later type than 
BLSl galaxies (Fig. [5] (left)). The two last plots illustrate 
again that NLSl bulges have less prominent (i.e. smaller, 
fainter, less cuspy) bulges than BLSl bulges. 

While the Sersic index is a convincing tool to distinguish 
pseudo-bulges from classical bulges ( Fishe r fc Dror y 2 00^ , 
iGadottil (|2009h uses the Kormendy relation {< fj.^ > —re) 
to identify pseudo-bulges as fainter bulges than predicted 
by the fundamental plane of elliptical galaxies. In Fig. [6l 
we present the Kormendy relation, the surface brightness 
magnitude at the effective radius, and the Sersic index versus 
the effective radius of the bulge. For Fig. [6] a) and b) we 
have overdrawn linear fits to the data for the BLSl sample. 
While these do not reveal any marked offset between the 
two NLSl and BLSl classes, NLSl bulges are systematically 
fainter than those of BLSl (i.e. they tend to lie under the 
line). This result is consistent with the common structural 
properties we are finding for the NLSl class. Finally, Fig. [6] 
c) shows NLSl and BLSl host bulges, together w ith pseudo- 
and classical bulges (from [Fisher fc Drory|[200^ ) in the rit, — 
re plane. Again, it shows clearly that NLSls lie in the region 
occupied by pseudo- bulges, while BLSl are spread over the 
whole range of pseudo- and classical bulge properties. 

Finally, we also look at the distribution of Sersic indices 
with the FWHM of the broad component of H/3 in Fig. (T] 
While the sample is not large enough to make conclusive re- 
marks, we note the existence of a correlation of the FWHM 
with the Sersic index, confirmed by a Spearman's rank cor- 
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Figure 6. Relations of bulge parameters with the effective ra- 
dius Te of the bulges, a) Kormendy relation, i.e. the mean surface 
brightness magnitude within the effective radius versus r^. b) Sur- 
face brightness magnitude at the effective radius versus r^. The 
solid lines are the linear fit found for BLSls from our sample. 
The magnitudes are given in the STMAG system, c) Effective ra- 
dius vs. Sersic index c omparing classical and pseudo-bulges from 
iFisher fc Drorvl ||2008| ) with our samples of NLSls and our BLSls. 
The NLSls seem to lie at the expected re — nt, of typical pseudo- 
bulges. 
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Figure 5. (Left) Distribution of bulge-to-total light ratio in NLSl and BLSls galaxies. According to fe.g. [Masters et al.ll20ich , NLSl 
galaxies would therefore be later type galaxies than BLSl. (Right) Verification of the mean B/T decrease with the Hubble Type, B/T 
against the bulge Sersic index and against the bulge effective radius. The two last plots confirm the link between B/T light ratios and 
the prominence of the bulge. 
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Figure 7. FWHM vs Sersic index. NLSls and BLSls are repre- 
sented by black and green dots respectively. The linear fit is based 
on the two samples. The Spearman's rank correlation coefficient 
is ~ 0.49. 

relation coefficient of ~ 0.49. This correlation indicates, at 
least in the low FWHM(iif/3) range, a possible connection 
between the bulge concentration and the broad line region. 

3.3.3 Complementary Studies 

Three other studies support our conclusion about the bulges 
of NLSls. 

In s tudying the central engines of NLSl galaxies, 
I Ryan et' al. (2007) perform a bulge disk decomposition of 11 
NLSls galaxies with z ^ 0.05 in the J-band and Ks-ha,nd us- 
ing adaptive optics data from the 3.6m CFHT. Their mean 
Sersic indices are < nj >= 1.52 and < tik >~ 1.38, and 
standard deviations (j„ of ~ 0.44 and ~ 0.48 respectively. 
As the Sersic index seems to be, at most, weakly corre lated 
to the photometric band (see e.g. iFisher fc Drorvll2008l . who 
compare th e Sersic index in the V and in the H bands), the 
iRvan et al.l decomposition supports our results^ 

In studying low black hole mass systems, iGreene et al.l 



(1200811 argue that most of their disk galaxies have pseudo- 
bulges. Since these systems have small FWHM(H;3), they are 
also likely to be NLSls although the relative Fell strengths 
they found are lower than in classical NLSls (| Greene fc Hd 
l2004h . 

In a more recent paper, iMathur et alj (|201ll ) study 10 
NLSl galaxies (ACS/HRC using the F625W filter). They 
perform a similar bulge disk decomposition (i.e. they fit a 
Sersic profile for the bulge and an exponential profile for the 
disk), and also conclude that they have pseudo- bulges. While 
our conclusions are based mainly on the S ersic index, the y 
use the Kormendy relation, as advised by iGad ottjl liooi), 
to conclude on the pseudo-bulge nature of NLSl bulges. In- 
deed, they do not find systematically low Sersic indices but 
obtain a mean Sersic index of < n >= 2.61 and a large Sersic 
index dispersion (t„ = 1.82. In fact only six out of the ten 
galaxies in their sample have Sersic index values consistent 
with pseudobulge profil es. To try and under stand this dif- 
ference, we note that the lMathur et af] l|201ll ') sample differs 
from that presented by us here. Specifically, their sample is 
found at larger redshift, < z >~ 0.24, leading to a scale 
in kpc/" on average ~4 times larger than in our sample. 
A large scale limits how well the fit is constrained by the 
central regions, and confusion in the light distributions be- 
tween the nucleus and the bulge may arise, possibly leading 
to higher Sersic indices. 

One caveat to these works is the lack of a comparison 
sample of BLSl hosts, which our results show is important. 
By including one, we show that one should consider the hosts 
of NLSls to be a sub-set of all BLSls, rather than being 
totally separate. The distinction between NLSl and BLSl 
hosts is thus that while NLSl hosts specifically have pseudo- 
bulges, BLSl hosts have a range of bulge types including 
pseudo-bulges, composite bulges, and classical bulges. 

3.4 Secular evolution has always prevailed 

The Sersic index distribution and the prominence of the 
bulge both indicate that, statistically, NLSl hosts have 
"pure" pseudo-bulges, in c ontrast to BLSl galaxies. The 
consequence of this result (|Kormendv fc Kennicutt] |2004| ) 
is that internal secular processes must have dominated the 
past evolution of NLSl hosts. And therefore it is from this 
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perspective that one should attempt to explain the particu- 
lar AGN properties observed in NLSl galaxies (such as low 
black hole mass, high accretion rates, etc.). 



4 NLSl EVOLUTION AND BLACK HOLE 
GROWTH 

We explore here the implications of our conclusion that sec- 
ular processes have dominated the evolution of NLSl galax- 
ies. We focus on the issue of the black hole growth of NLSls, 
and in understanding whether or not NLSl galaxies are in 
a special phase of black hole growth. 



4.1 Expected galaxy populations that have 
evolved through secular evolution 

Over the last few years, the relative importance of major 
mergers versus minor mergers and secular processes in driv- 
ing galaxy formation and evolution has become a key issue 
in s imulations and sem i-analytic models. iGenel et all (|2008l ) 
an d IParrv et all (|2009l ) offer two different perspectives to un- 
derstanding the growth processes of dark matter halos and 
galaxies. Both studies conclude that major mergers are not 
nec essarily the main dr iver of galaxy mass evolution. 

iGenel et al. I (|2008) investigate, by analysing cosmolog- 
ical simulations, the growth of dark matter halos. They ex- 
tract halo merger fractions and mass accretion rates from 
the Millenium simulation in order to study the possible role 
of major mergers in the evolution of halos from z ^ 2 
to z — 0. Following the fate of halos in the mass range 
11.5 < logM^=2.2 sS 12.8, they find that ~l/3 of halos 
which reach z—Q have not undergone any major mergers 
since z ~ 2.2 and that such halos gain > 70% of their new 
mass via mergers less intense than 1:10, demonstrating the 
import ance of non - major merger processes. In a following 
paper, iGenel et al.l (|2010l ) also show that, independently of 
halo mass, ~ 40% of the mass in halos have been assembled 
throug h smooth accret ion. 

Parrv et al.l (|2009l ) study galaxy growth by analysing 
two different gala:xy formation models, both also based on 
the Millenium simulation. Their statistical results are reveal- 
ing. For both models, they find that only ^49% of ellipti- 
cals, ^3% of SOs and ^2% of spirals undergo a main branch 
major merger (mass ratio greater than 1:3) in their entire 
formation history. In other words, ~98% of spiral galaxies 
- which are the most common morphological type of NLSls 
host galaxies - do not undergo any major merger. These 
results are largely independent of total stellar mass of the 
galaxy except for ellipticals. In a further step, they quantify 
the relative impact of disk instabilities, major mergers and 
minor mergers on galaxy morphology by determining the 
stellar mass fraction from each process as a function of the 
total stellar mass. For both models, they find that instabili- 
ties and minor mergers are the main mass contributor, with 
their relative contributions dependin g on the treatmen t of 
disk instabilities in the models (see Parrv et al.ll2009l and 
references therein). 

While it is beyond the scope of this paper to discuss the 
treatment of the various physical processes in these galaxy 
formation models, one clear conclusion is that even hierar- 
chical cosmological simulations give rise to large galaxy pop- 



ulations that have evolved through secular processes. Inter- 
esti ngly, observational stu dies also reach similar conclusions. 

IWeinzirl et~al] l|2009l ') study present day (D < 60Mpc) 
spiral galaxies. By performing two dimensional multi- 
component decompositions of 143 high mass (M ^ 10^° M0) 
spirals, they analyze the bulge Sersic index and B/T distri- 
butions. Their results highlight the large fraction of bright 
spirals having B/T<0.2 (~69%) and n ^ 2 (~76%) where 
many of them host bars (~66%). By comparing their result 
to theoretical predictions, they find that ~ 66% of present 
day high-mass spirals have not undergone a major merger 
since z ^2 and likely not even since 2^4. This conclusion 
conveys the importance of minor mergers (in the present 
case for a mass ratio < 1 : 4) and secular processes since 
z < 4. 

Finally, 'Cister nas et^ (|2010l ') recently analyzed the 
relevance of different triggering mechanisms for AGN activ- 
ity. Based on visual analysis of 140 AGN and 1264 inactive 
galaxies with HST imaging, they measure the fraction of 
distorted morphologies which they take to be a signature of 
recent mergers. They conclude that the bulk of black hole 
accretion has been triggered by secular processes and minor 
interactions since z ^ 1. 

By assessing the role of secular processes, these various 
theoretical and observational lines of evidence offer a cosmo- 
logical context to our conclusions concerning NLSl hosts. 
They show that a significant fraction of galaxy mass and 
a large number of galaxies have evolved from early cosmic 
times without any mergers. It is reasonable to suggest that 
NLSls, which our analysis shows must have evolved without 
mergers, may be tracers of this population of galaxies. Our 
attempt, in the rest of this section, to put this hypothesis 
on a more quantitative basis leads us to an estimate of the 
duty cycle of NLSls. 



4.2 How common are NLSls? 

In the last decade, various surveys of nearby galactic nuclei 
have found that the fraction of objects classified as AGNs 
is surprisingly large. They show that, of all local galaxies, 
approxima tely 10% are Seyferts and 40% can be considered 
active (see iHolboOSl and references therein). 

Looking at the NLSls, several surveys using opti- 



cal and X-ray selected samples (e.g 



ICrenshaw et al.ll2003l : IZhou et al]|2006l ') find that thev make 



Williams et all I2OO2I : 



up approximately 15% of Seyfert 1 galaxies. Based on the 
unified AGN scenario, one can expect that this fraction 
should apply also to type 2 Seyferts: that ~15% of Seyfert 
2 galaxies may have narrow broad emission lines that are 
hidden from sight by the obsc uring torus. Such a pos sibilty 
has already been suggested bv lZhang fc Wang (|2006l ). who 
argue that Seyfert 2s without a hidden BLR (i.e. one that 
cannot be observed in polarised light) are the counter-parts 
of NLSls. Taking this further, one could reasonably argue 
that the NLSl definition might extend even to the lower 
luminosity AGN. 

Combining these fractions together leads to the assess- 
ment that 2-6% of local galaxies could be "NLSl-like" galax- 
ies, i.e. active galaxies with "narrow" broad lines whether 
obscured or not, the evolution of which has been dominated 
at all cosmic times by secular processes. 
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4.3 Duty cycle of NLSl black holes 

Our analysis in section [3] shows that NLSl galaxies sta- 
tistically have a bulge Sersic index nt < 2 and a bulge- 
to-to tal light ratio B/T < 2. The observational study of 
IWeinz i rl ct a l. (2009), corroborated by the theoretical study 
of iParrv et al] (|2009l '). argues that approximately 2/3 of lo- 
cal spiral galaxie^ have similar properties to NLSl hosts. 
Assuming that each galaxy among this population can po- 
tentially undergo nuclear activity and become "NLSl-like'jf|, 
we can estimate the duty cycle of NLSls. We argue above 
that "NLSl-like" objects may comprise as much as of 
local galaxies (note that we use the upper end of the range 
above in order to be conservative later in our estimates of 
black hole growth). This implies that their duty cycle should 
be around ^^9%. 

Since the Hubble se quence formed at z^l 



J van de n Bergh et all |2000|: iKaiisawa fc Yamadal l200ll: 
Conscli ce et al.l |2004| : lOesch et al.l l2010l : ICassata et"aLl 
2010l) with bars be c oming numerous at t his redshift 
jAbraham et al. 19991: Elmegreen et al.l 20041 : Jogee et al.l 



12004 iMarinova fc Joged l2007l : ISheth et al.l 120031 , 120081) , 

accretion driven by large-scale bars can only have occurred 
in the last ~7.7Gyr. With a duty cycle of ~9%, this means 
that NLSls have actively accreted onto their BHs for 
~690Myr. 

Assuming their black holes are accreting at the Edding- 
ton rate, the e-folding time of their BH build-up is given by 
the Eddington time-scale, ts ~ 4.4 x 10^ yr. Therefore, the 
black hole mass increase is given by (e.g. IVolonterillioia i : 



MBH(t) 



M^'!^'^{to) exp f ^ 



tE 



(1) 



where t is the current (observed) time, and t — t ~ to is 
the total accretion duration since the initial time to- The 
radiative efficiency e is a key parameter, and can have a 
major impact on the black hole growth rate because it is 
inside the exponential term. The standard value of e is ~ 0.1. 
Adopting this then, and assuming a seed mass A{§'^'^{to) — 
10'' — lO'* Mq, the current mass AfBH(t) could range from 
10^ to 10^°Mq. On the other hand, accretion onto a fast 
rotating Kerr BH can lead to a radiative effi ciency e ~ 0.2 
or higher l|Volonteri et al.ll2005l : IJogee II2006I ). For the same 
M^^'^{to), this higher radiative efficiency implies a much 
slower BH growth and leads to current BH masses between 
5 X lO'' and 5 x lO'^ Mq. 

Based on these very simple estimates, we reach two con- 
clusions. Firstly, NLSls are not necessarily in a special phase 
of their black hole growth. Their black holes have required 
7-8 Gyr to grow to their current size. To increase their black 
hole mass by another factor 10 requires - for a 9% duty cy- 
cle and e — 0.2 - another ~2.8 Gyr. Thus, despite their high 
Eddington ratios, NLSls are not imminently evolving into 
BLSls, although they may do this eventually. Secondly, the 
low BH masses of NLSls (with respect to BLSls) can easily 



^ high mass (M* > 1.0 X 10-'^'' Mq) low-to-moderately inclined 
(i < 70°) spirals. 

clearly the 2/3 of the local galaxies are not necessarily potential 
true NLSls; but, given that their hosts appear to have evolved 
over cosmic time in a similar way, they might be "NLSl-like", as 
defined in Section 14.21 during their accretion phase. 



be accounted for by a high radiative efficiency. The theoret- 
ical and observational evidence for rapidly spinning black 
holes in NLSls is the topic of the section [STTI 



5 DISCUSSION 

5.1 NLSls and highly spinning black holes 

As discussed in the last paragraph, the radiative efficiency 
of the AGNs have a direct impact on their black hole growth 
and their final black hole masses. But radiative efficiency is 
determined by the supermassive black hole (SMBH) spin, 
which in turn is influenced by the AGN accretion his- 
tory. Different scenarios have been studied where SMBH 
spins evolve through mergers, subsequent prolonged accre- 
tion (const ant angular momentum axis of the accreting m a- 
terial, e.g. IVolonteri et alllioosl : iBerti fc Volonterill2008l l or 
chaotic accretion (rand om disc orientations with respect 
to the black hole, e.g. iKing fc Pringld [20071: iKing et all 
20081: iBerti fc Volontcri 20081). IVolonteri et al.l (2005) and 



Berti fc Volonteri (2008) have shown that if major and mi- 



nor BH mergers are the sole source of material, then the 
distribution of BH spins in a z = galaxy population will 
reflect that of the initial seed BH spin. In contrast, they find 
that prolonged gas accretion triggered by galaxy mergers 
tends to spin BHs up, and that galaxies where a significant 
fraction of the BH growth occurs in this mode could have 
maximal BH spin. However, in such a case, most distant 
quasars would have high radiative efficiency and would in- 
efficiently grow their black hole. This would require massive 
BH seeds, in conflict with the Soltan argument. If instead, 
the accretion proceeds by short rando mly oriented events 
(jKing fc Pringldl2007l : [King et al.ll2008l ). then the spins wiU 
tend to be low and lead to high BH masses, resolving the 
conflict with the Soltan argument. But NLSl BHs are likely 
fed via secular rather than merger processes, therefore the 
angular momentum of the infalling matter could be related 
to the host structure and hence have a favoured direction. 
NLSl SMBHs could be evolving through the prolonged gas 
accretion scenario. Therefore, NLSl secular evolution would 
imply high spins and low BH masses. 

From an observational point of view, it appears to be 
possible to derive BH spins using the Fe Kq line in the 
hard X-ray continuum using accretion disk reflection mod- 
els. While there are several potential reasons that might pre- 
vent one from obtaining a good spin constraint for a given 
source (such as too few photons, or too narrow iron line 
hamperin g the flt to pick out the role of relativistic con- 
tributors, |BremMman[[2003) several recent works have been 
able to derive formal const raints on BH spin for a num- 
ber of Seyfert galaxies (e.g. 'Brenneman"2007'; 'Fabia n et al.l 
[2009; Miniutti ot al. 2009; Ponti ct al. 2010). Among these 
measurements, we find that most of them are for NLSls, 
and that the derived spin values are very high. While there 
are still too few measurements to draw any firm conclusion, 
these results are already suggestive that NLSls could have 
highly spinning BHs owing to the prolonged disk accretion 
onto their BHs. 

Finally, as we have already touched on in section a 
high spin leads to a high mass-to-energy conversion or radia- 
tive efficiency (because the last stable orbit is closer to the 
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horizon of the BH), and hence to a slow BH growth. There- 
fore, the high Eddington ratio typicaUy observed in NLSls 
could be a signature not just of highly accreting BHs, but 
also of rapidly spinning BHs. Given that pseudo-bulges are 
one consequence of secular evolution in a galaxy, and that 
another consequence is that the SMBH should be rapidly 
spinning, our result that NLSl hosts have pseudo-bulges 
leads to the prediction that a large fraction of NLSl BHs 
should have very high spin. Future X-ray missions may be 
able to test this. 



5.2 Black hole - bulge scaling relation 

One question concerning NLSls that has received much at- 
tention is whether they follow the Mbh — cr* relation or are 
offset under it. The Mbh — cr*, or more generally the BH- 
bulge scaling relations, are often interpreted as physical evi- 
dence for the co-evolution of the central BHs with the galac- 
tic bulges. The case of NLSls is rather controversial. On one 
hand, many studies suggest they may reside below the re- 
lation, in which ca se they could be evo l ving onto it (e.g. 
Mathur et all I2OOII: iBi an & Zhao ' 2004 iGrupe fc Mathuil 



2004 iMathur fc Grupeli2005a .bi: iZhou et all 120061 ). On the 



other hand, several different studies place them on the rela- 
tion once contaminating effects have been correc ted for, such 
as [O III] l injfl broadenin g due to outflows (see iBotte et al.l 
I2OO4I andU omossa fc X u 2007), or radiation pressure (as 
proposed bv lMarconi et al.i,200S ). 

Current developments regarding the Mbh — cr* rela- 
tion highlight that it may not be universal (common to 
all morphological types), and that perhaps one should dis- 
tinguish betw een barred and barless galaxies, disks and el- 
lipticals (e.g. iGrahamlliooi IHuI|2008I : iGraham fc Lill2009l : 



iGraham et al 



201 ih. c lassical bulges and pseudo-bulges 



(e.g. iNowak et al.ll201ol : [s ani et al. l l2011xl : lKormendv et alj 
[201^. In fact, barred, disk galaxy bulges and pseudo-bulges 
appear either to lie below the relation or to scatter around 
it. Additionally, on a more theoreti cal side, some authors 
l|Penall2007l : Ijahnke fc Macciol [ioill ) suggest that the BH- 
bulge scaling relations could be non-causal (their origin 
would not invoke a physical coupling between the SMBH 
and the galaxy) but rather would be naturally produced by 
the merger-driven assembly of bulge and BH masses, and 
therefore galaxies with pse udo-bulges would not be expected 
to obey the same relation ijjahnke fc Macci6ll201ll ). 

Our result that NLSls have pseudo-bulges suggests that 
we should expect these AGN to lie in the same region as in- 
active galaxies with pseudo-bulges, that is scattered around 
and below the Mbh — cr* relation. It is not yet understood 
how - or whether - black hole and bulge growth are linked 
when secular processes drive their evolution. Thus, while it 
is clear that their black holes are still growing, we cannot 
predict where they will end up on the Mbh — o"* plane. 



5.3 Evolutionary scenarios 

Several authors have suggested different links between NLSl 
galaxies and other AGN types in evolutionary sequence con- 
texts. We briefly discuss them in the light of our results. 

iMathuil (|200G| ) argues that NLSl galaxies might be in 
an early stage of evolution owing to their small growing 
black holes and higher Eddington rates. This proposition is 
not inconsistent with our results. Nevertheless, it illustrates 
a different perspective. Either NLSls would have their nu- 
clear supermassive black holes recently form ed and NLSls 
would be young ob jects evolving into BLSls l|Mathuill2000l : 
iMathur et a"i]|201lD . or NLSls would not be in any special 
phase of their evolution but simply have BHs that are grow- 
ing slowly due to their duty cycle and spin. However, both 
perspectives agree that NLSls galaxies have pseudo-bulges 
and that their black hole growth is driven by secular pro- 
cesses as opposed to mergers at high redshift. 

iKawakatu et al.l |2003) proposed an evolutionary track 
from type 1 ULIRG, to NLSl to BLSl. The connection be- 
tween ULIRG and NLSl appears contradictory with our re- 
sults. Indeed, local ULIRGs are the result of galaxy merg- 
ers, while we have argued that, based on the host properties, 
NLSls have a secular driven evolution. 



^ this line is often u sed as a surrogate of the stellar dispersion 
Oi,, see lNelsorJ 1I2OOOI ). 



IZhang fc Wand (|2006| ) study Seyfert 2s with and with- 
out a hidden BLR (i.e. presence or absence of BLR in polar- 
ized light) and suggest that non-HBLR Seyfert 2s are the 
counterparts of NLSls vie wed at high inc l inatio n angles. 
In their subsequent paper, IWang fc ZhangI (|2007t ) propose 
an evolutionary sequence of the narrow objects to broad 
line AGN considering time evolution of the black hole mass 
and the accretion rates. While it is not inconsistent with 
our results, the distribution of bulge properties in BLSls 
suggests that not all BLSls come from NLSls that have 
evolved secularly, but that the BLSl population should in- 
clude galaxies that have undergone interactions and mergers. 
While IWang fc Zhan3 f2007. ) propose a secular evolution 
from NLSls to BLSls (NLSls would be an early AGN phase 
and would evolve to BLSls during the AGN activity time), 
IZhu et ai] (|2009l ) proposed a similar scenario, but where 
NLSls would be produced by mergers of smaller galaxies 
compared to BLSls and could evolve to BLSls only if they 
encounter more mergers to grow them. This last scenario 
appears contradictory to the results of the present paper. 

While at this point there is no consensus on the cos- 
mic evolution of NLSls, our results suggest that they are 
a special case in which the evolution has been dominated 
at all time by secular processes. Thus, if the BHs in NLSls 
continue to grow, they must eventually become broad-line 
AGN, and thus become part of the BLSl population. How- 
ever, our results also show very clearly that not all BLSls 
have grown in this way, and that mergers have played a role 
in the evolution of the BLSl population. In this respect, per- 
haps the most enlightening question would be: when NLSls 
evolve into BLSls, will they be distinguishable from systems 
classified as BLSls but having undergone galaxy interactions 
and mergers? Perhaps one can already begin to address this 
by studying the BLSls with pure pseudo-bulges, and asking 
whether they have definable characteristics that differ from 
the BLSl population as a whole. 
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6 CONCLUSIONS 

From a review of the literature, we show that secular evolu- 
tion in NLSl galaxies is a powerful and on-going process on 
all scales, in contrast to BLSl galaxies. To assess the role 
of secular processes in the past evolution of NLSl galaxies, 
we examine their bulge properties by performing bulge-disk 
decompositions on NLSl and BLSl galaxies with archival 
HST images. The results indicate that NLSl host bulges are 
pseudo-bulges and distinct from the much broader popula- 
tion of BLSl bulges. Prom these results, we conclude that 
NLSls represent a class of AGN in which the black hole 
growth is, and has always been, dominated by secular evo- 
lution. 

Such am evolutionary mode signifies also a different 
black hole growth mode in NLSls than in merger-built sys- 
tems. Interestingly, simulations of prolonged disk-mode gas 
accretion onto black holes show that the most efficient way 
to spin-up a black hole is through smooth accretion of mar 
terial. In this light, our results suggest that NLSl galaxies 
should possess highly spinning BHs which is indeed what 
has so far been observed. 

Our picture of the NLSl galaxy phenomenon can be 
expressed as follows. The activity in NLSl galaxies is, and 
always has been, powered by internal secular processes. This 
has lead to the growth of a pseudo-bulge. It is also character- 
ized by a disk-mode accretion onto the central object, which 
tends to spin up the black hole. This leads to high radiative 
efficiency of the accreting material, therefore reducing the 
actual mass accreted onto the black hole and slowing its 
growth. The high radiative efficiency could in part explain 
the high Eddington ratios and small black hole masses of 
NLSls. 

This picture can be tested by analyzing the angular mo- 
mentum in NLSl bulges to assess definitively their pseudo- 
bulge nature. And studying the kinematics in the central 
part of NLSls would help to understand how gas is trans- 
ported to their central regions, what the mass inflow rates 
are, and the role played by angular momentum. In parallel, 
systematic measurements of black hole spins by future X-ray 
missions would also shed light on the growth of their black 
holes at the smallest scales. 
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APPENDIX A: SERSIC LIGHT PROFILE 

The Sersic power-law intensity profile is frequently used in 
the study of gal axy morphology. I t has the followi ng func- 
tional form (e.g. ICaon et al.ll 19931 : IPeng et al.ll2010l ') 



/(r) = le exp 



l/n 



(Al) 



where 7^ is the surface brightness at the effective radius r^. 
The parameter re is known as the effective, or half-light, 
radius, defined such that half of the total flux lies within Ve. 
The parameter n is the Sersic index: with n = 4 the profile 
is the de Vaucouleurs profile typical of elliptical galaxies, 
n — 0.5 gives a Gaussian, and n = 1 is the exponential 
profile typical of disks. Finally, the parameter k is in fact 
coupled to n and is not a free parameter. 

The exponential disk profile is more frequently de- 
scribed by the compact form 



I{r) = Jo exp ( -— ) , 



(A2) 



where Jq is the central surface brightness, Jq = /efi", and 
the scale radius given by re — 1.678ri,, or more generally by 

Te = k'Vs. 



Al Sersic luminosity ratios 

The bulge-to-disk (B/D) and bulge-to-total (B/T) luminos- 
ity ratios are relevant quantities for the study of galaxy 
morphology. These ra tios tend to d e crease from early- to 
late- type spirals (e.g. [Masters et al.1 (jloi^))- The B/D ra- 
tio, where the disk is described by an exponential profile, 
can be expressed analytically as 



B ^ n,r(2n6)e" f qtRl\ fh 
and the B/T as 



B ^ ntr{2nt)eyti^''»qtRlle 

T ni,r(2ni,)e«/K2"i.g6i?2/, + qaRpo ' 



(A3) 



(A4) 



where the subscripts b and d refer to the bulge and the disk 
respectively. le is the effective surface brightness of the bulge 
and Jo is the central surface brightness of the disk, Ri, is the 
effective radius of the bulge and Rd the scale radius of the 
disk, and qt, qd are the axis ratios of the respective profiles. 
To calculate the p arameter k, we use the analytic expansion 
Equation 18 from lCiotti & BertinI (|l999l ) valid for n > 0.36. 
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Step 1: Fit of 
the Disk 



Step 2: Fit of , 

the Disk+Bulge 

i 

Step 3: Fit of ^ 

the Disk+Bulge+Nucleus 



Step 4: 

Iteration and Selection 



Initial guesses 

for disk component. 



Fit of the exp. profile. 
Central region masked. 



Initial guesses 
for the bulge. 
Output of Step 1 
for the disk. 



Fit of an exp. and 

Sersic profiles. 

Small core region masked. 



Jleraiion unlil 
all param. 



Initial guesses 
for the nucleus. 
Output of Step 2 for 
the disk and the bulge, 



Fit of an exp., a Sersic and 
a gaussian or PSF profiles. 



] Iteration until 
all param. 
free to vaiy 



Choose the best fit from Step 3 based on: 

(1) (radial) chi2, (2) (radial) residual, (3) radial profile 
Parameters to vary: 

(1) regions masked, (2) sky value, (3) psf 



Figure Bl. Our iterative fitting procedure. 

APPENDIX B: FIT PROCEDURE 

To accomplish the bulge-disk decomposition, we use the two- 
dime nsional profile fitting algorithm Galfit (jPeng et alJ 
which aUows, by minimizing the value, to model 
the light profiles using e.g. Sersic profiles. 

Since we fit HST images, we generate the PSFs with 
the TinyTiwQ code for the WFPC2 camera as well as the 
F606W filter. For each image, we create a PSF referenced 
to the galaxy center with an uniform weight along the wave- 
length range. 

Initial parameters are estimated from visual inspection 
of the images (positions, P. A., ellipticities, radii, etc.). Nev- 
ertheless, our iterative fit procedure, as described hereafter, 
ensure that the choice of initial estimates do not influence 
the final results. 

Our procedure can be divided in 4 steps as summarized 
in Fig. IBll For each step we create appropriated masks to 
remove regions from the fits. 

1. We start by fitting only the disk of the galaxies using 
an exponential profile and masking the central region. 

2. Once a reasonable model is obtained, we fit an ad- 
ditional Sersic component to model the bulge. Our initial 
estimate of the Sersic index is nt, = 2 (i.e. the index thresh- 
old between pseudo- and classical bulges), while the other 
parameter initial estimates are based upon visual inspection. 
When fitting, we first keep the outputs of our step 1 fixed in 
the fit and then free to vary (together with the Sersic profile 
parameters). During this step we mask only the core of the 
galaxies (5 to 10 pixels radius) to avoid possible influence of 
the central nucleus or AGN source. 

3. We then fit an additional Gaussian (with initial 
FWHM of 2.5 pixels) to model the AGN and iterate, if 
necessary, until all parameters are free to vary. As a non- 
negligible fraction of the images presents a saturated core 
with charges leakage along the columns (34% in total, 50% 
in the NLSl and 26% in the BLSl samples), we mask care- 
fully these pixels. We also verify the resulting FWHM of the 
Gaussian and its axis ratio. An additional motivation to use 
a Gaussian instead of a PSF is to model any nuclear star 
cluster which, if not accounted for in the fit, will artificially 
increase the bulge Sersic index. 

During these three steps, we judge the quality of the fit 
based on the residual images, the values, and the param- 
eters values (ensuring that they are physically meaningful : 
the magnitude of the components have to be greater than 

^ http:/ /www. stsci.edu/software/tinytim/ 



the sky level, the Sersic index value should be acceptable, in 
particular nt < 8, and the physical radii -Ri, and Rd should 
be reasonable). 

4. At the end of the third step, we additionally look 
at the radial profiles of the model, its components, and the 
original image. These plots, together with the radial residual 
and the radial enable us to diagnostic possible problems 
and identify influences of non-fitted structures. Upon exam- 
ination of these plots and the image residuals, we decide 
whether it is necessary to mask relevant structures such as 
rings or spirals. We refine the masks and the fits until the 
radial profile of the model does not appear to be influenced 
by these structures but do translate correctly the bulge and 
the disk components. 

While the Galfit outputs provide errors, those are 
purely statistical and are insignificant compared to system- 
atic errors. In the following sections, we discuss the robust- 
ness of our fits to different parameters: the treatment of 
additional structures (deviation from idealized profiles), the 
unprecise knowledge of the background level, the core sat- 
uration and the PSF. Those parameters give indirect infor- 
mation on the systematic errors. Finally, we give in Fig. IB2I 
and Fig. IB3I eight examples of our fits where we indicate by 
open symbols the radial ranges entirely or partially masked. 
The radial residual and the radial xt, given in the upper 
panels, also provide indirect information on the fit errors. 



Bl Choice of fitting range: treatment of 
additional structures 

As already mentioned, any nuclei is accounted for in the 
fits by the use of a small Gaussian profile. The alternative 
consists of removing the nucleus from the fit by masking 
it. Nevertheless, with such a procedure, the bulge-disk de- 
composition is sensitive to the quality of the mask of the 
nucleus and can also be affected by the reduced number of 
constraints, i.e. the bulge can be underconstrained if the 
mask is too large. Despite these considerations, we try the 
alternative and mask systematically the central region of 
the images. We use a default circular mask with 10 pix- 
els radius, adjusting it only in a few cases (4/28) to have 
a minimum radius of 5 times the standard deviation o of 
the Gaussian component (2/4) or to keep a reasonable num- 
ber of constraints for the bulges (2/4). The median Sersic 
index difference between the fits with nucleus masked and 
the fits with the nucleus included is reasonably small: 0.22 
and 0.29 for the NLSl and BLSl sample respectively. This 
test shows that to mask the nucleus instead of fltting it by 
a Gaussian component biases slightly the Sersic indices to 
higher values. Nevertheless, except for a few cases (3/28 with 
Arifc > 2 and 2/28 with 2 > An;, > 1), the increase remains 
small, acknowledging the robutness of the flts. 

As our samples are made of high resolution HST data, 
they also present many structural details. The effect of de- 
viations of surface brightness such as rings, bars or spirals 
are considered in the last step of our fitting procedure (step 
4). Upon examination of the radial profiles and the resid- 
ual image, we identify, if any, potential additional struc- 
tures and manually create appropriate 2D masks. We then 
refit our model and iteratively adjust the mask according 
to the radial profile and residual image. We recognize that 
this practice is subjective but nonetheless we believe it to 
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Figure B2. Bulge-Disk decomposition illustration from the NLSl sample. The radial ranges of the mask are indicated by open symbols 
and straight lines in the radial profiles. The upper plot gives the radial distribution of the reduced x^i the middle plot is the magnitude 
difference Afi between the original and the model images, the bottom plot presents the observed (dots), the total model (solid line), the 
modelled bulge and the modelled disk (dashed) light distributions. The radius is given in r^/** to emphasize the central region. 



be necessary. Indeed, modelling structural details - such as 
inner rings - is a complex task, and the fit of any addi- 
tional component requires more constraints, which cannot 
be provided by our single snapshot HST images (indeed the 
current bulge/disk models already give ~ 1 for most of 
the fits). The alternative practice of not masking these struc- 
tures could lead to wrong results. For example, in the case 
of MRK 42, presented in Fig. IB2I we first do not mask the 
inner-ring and obtain a compact bulge with nt ~ 0.7 (in our 
step 3), while after masking the stucture, we obtain a more 
reasonable Sersic index nt ~ 1.27 (in our step 4), which - 



according to the radial profile in Fig. IB2I - is a much more 
accurate model of the bulge. Finally, our iterative process 
ensures that the regions we mask are physical additional 
structures in the galaxies. 

B2 Background level 

Given that the background may be coupled to the disk pro- 
file, an important aspect of bulge-disk decomposition is to 
correctly fix the sky level. Since no precise sky measurements 
are available for our HST images and that the galaxies filled 
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Figure B3. Bulge-Disk decomposition illustration from the BLSl sample. The radial ranges of the mask are indicated by open symbols 
and straight lines in the radial profiles. The upper plot gives the radial distribution of the reduced x^i the middle plot is the magnitude 
difference A/x between the original and the model images, the bottom plot presents the observed (dots), the total model (solid line), the 
modelled bulge and the modelled disk (dashed) light distributions. The radius is given in r^/** to emphasize the central region. 



most of the field of view, we have to fit the background 
level together with the other components and minimize the 
couphng with the disk profile. Therefore, we take a partic- 
ular care to any additional light profiles in the images such 
as stars, satellite galaxies or remaining large cosmic rays. 
These are carefully masked so to minimize their influence 
on the sky level. Despite this particular attention, we fix the 
sky level in 2 cases - MRK279 and MRK609, both part of 
the BLSf sample - at the value obtained in our step f (i.e. 
fit of the disk and the background only). If the sky is not 
fixed, the resulting parameters are not physical: in the case 



of MRK 279 the background becomes negative and the ra- 
dius of the disk excessively large; in the case of MRK 609 
the sky becomes extremely large and the disk and the bulge 
shrink (with nj, < 0.5). 



APPENDIX C: RELIABILITY OF THE 
BULGE/DISK DECOMPOSITION 

We analyze here possible effects influencing our fit results 
the core saturation in the images and the choice of PSF. 
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As already mentioned, a non- negligible fraction of the 
images presents saturated core with charges leakage: 50% in 
the NLSl sample and 26% in the BLSl sample. To minimize 
any effect on our fits, we mask carefully the saturated pixels 
and the pixels affected by charges leakage. Looking at the 
Sersic indices, we find that gala:xies with saturated core have 
a mean Sersic index < ni, >= 2.27 and galaxies without sat- 
urated core have < rib >= 2.06. Therefore, if the saturation 
still affect our images, it would tend to increase the Sersic 
index. As our NLSl sample is more affected by saturation 
than our BLSl sample, it would tend to increase the mean 
Sersic index of NLSls. Consequently, any remaining effect 
of the saturation cannot account for the difference between 
NLSl and BLSl host bulge Sersic indices but would tend to 
decrease it. 

We test the dependence of our fit to the PSF used 
for convolution by refitting our NLSl sample with differ- 
ent PSFs. The PSFs for this test arc also generated with 
TinyTim but, instead of using an uniform weight along the 
wavelength range, they arc produced at the central mono- 
wavelength of the filter F606W. The difference in Sersic in- 
dex ranges from to 0.01 except for one object where it is 
0.06. Therefore, we conclude that the PSFs arc not critical in 
our fits and that our choice of using an uniform weight along 
the wavelength range to create the PSFs with TinyTim is 
acceptable. 



